Correlation lengths for vortex-liquid freezing 
in a model of layered high-temperature superconductor 



A. Zamora* 

Institute de Ciencias Nucleates, UNAM, Apdo. Postal 70-543, Mexico D.F. 04510 

(Dated: March 4, 2009) 

We present results from extensive numerical simulations on the Lawrence-Doniach model within 
the lowest Landau level approximation in the unconventional spherical geometry. We study the 
in-layer pancake vortex density- density correlation function (intra-layer structure factor) and the 
layer-to-layer order-parameter correlation function along the direction perpendicular to the layers. 
Our results show strong evidence for the existence of a single first-order phase transition at which 
both inter-layer coupling and appearance of crystalline order in the in-layer vortex correlations take 
place as temperature is lowered. 
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The study of the phase diagram of layered super- 
conductors in the presence of a magnetic field has at- 
tracted much attention since the discovery of the high- 
temperature superconducting cuprates 1 . In particular, 
the most interesting effects are found in the mixed phase 2 
where the external magnetic field penetrates the sample 
in the form of quantized flux- lines (or vortices). Within 
the mixed phase, originally described by Abrikosov as 
a vortex lattice, a high-temperature phase with finite 
resistivity has been found and called a vortex liquid. 
While the vortex-liquid phase is quite narrow in con- 
ventional type-II superconductors, its width is known 
to become large in highly anisotropic cuprates 4 (e.g. 
Bi 2 Sr 2 CaCu208.) Assuming a specimen completely free 
of defects, the low-temperature phase is expected to be 
the Abrikosov vortex-lattice and the phase boundary sep- 
arating the vortex liquid and vortex lattice appears to be 
a line of first-order phase transitions usually called the 
melting line 5 (which essentially coincides with the mag- 
netic irreversibility- line 6 ). Results from experiments, an- 
alytic calculations, and numerical simulations on several 
models of layered superconductors have shown that, in 
addition or simultaneous to melting, there is a transition 
associated with layer decoupling 7-13 . In another simu- 
lation, Wilkin and Jensen 14 on a layered model find a 
first-order transition associated with layer decoupling fol- 
lowed by melting but which shows no obvious thermody- 
namic signature. An alternative scenario was proposed 
by Kienappel and Moore 15 . From Langevin-dynamics 
simulations on the Lawrence-Doniach (LD) model within 
the lowest Landau level (LLL) approximation and in the 
same spherical (S) geometry that we consider here (we 
shall call this the S-LD-LLL model) they find hysteretic 
effects on several quantities which disappear on increas- 
ing the interlayer coupling-strength. They attribute this 
result to a line of first-order transitions ending at a crit- 
ical point which separates phases of coupled and decou- 
pled vortex-liquid (expected to freeze only at zero tem- 
perature.) 

In this paper we present an alternative view based on 
Monte Carlo (MC) simulations in the same S-LD-LLL 
model. From studies of inlayer vortex-density correla- 



tions as well as layer-to-layer order parameter correla- 
tions we find strong evidence for vortex-liquid freezing 
(simultaneous to inter-layer coupling) at a finite temper- 
ature. As in Ref. 15 , we model this layered superconductor 
by a system of M concentric spherical-layers of thickness 
do *C R, where R is the average radius of the sample, and 
interlayer spacing s > cIq. The reason for the choice of 
this geometry, instead of the usual array of plane-layers, 
is that this geometry guarantees full rotational and trans- 
lational invariance of the vortex system (which quasi pe- 
riodic boundary conditions in the planar geometry do 
not). This allows for the in-layer free movement of vor- 
tices without any boundary constraint, which is expected 
to better describe a real superconducting sample. This 
has been previously discussed in detail for the single-layer 
case by Dodgson and Moore 16 . An external magnetic 
field H(r) = (H R 2 /r 2 )r [with H a constant field], or- 
thogonal to all layers in every point, is generated by an 
infinitely long and thin solenoid ending at the center of 
the concentric spheres. We neglect the fluctuations in the 
magnetic induction inside the sample so that it has the 
mean value B = ^H{R) across all superconducting lay- 
ers. Flux quantization requires that the magnetic induc- 
tion penetrating the sample be such that 4irR 2 B = NQq, 
where $o = h/2\e\ is the flux quantum and N the num- 
ber of flux lines or vortices. Our choice of gauge that 
satisfies B = V x A is A = (n H R 2 /r) tan(6»/2)</>. We 
assume that the system without impurities accepts a de- 
scription in terms of the LD Hamiltonian functional 17 
which, under the previous assumptions, can be written 
as 
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Here ip n ( r ) is the two-dimensional Ginzburg-Landau 
(GL) order parameter in layer n and = — ih^ ± — 2eA 
the gauge covariant derivative operator (acting on the 
surface of the sphere). a(T) and (3 are the usual param- 
eters from the GL theory 18 . m± and mil are the effective 
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Cooper-pair masses in the perpendicular and parallel di- 
rections relative to the radial magnetic field. The LLL 
approximation consists in expanding the GL order pa- 
rameter in each layer as a linear combination of the de- 
generate eigenfunctions of D\ = Dj_ • Dj_ corresponding 
to the lowest eigenvalue 2|e|?i^o-ff. An orthonormal set 
of functions for the LLL on the surface of the sphere is 19 

<p q , N (0, <f>) = N q , N e^ sin"(0/2) cos JV -«(^/2), (2) 

where q = 0, 1, . . . , N labels the degeneracy of the LLL. 
The normalization factor is A/" r>s = [(s + l)CJ/47ri? 2 ] 1/2 , 
with CJ = s!/H(s — r)! the binomial cocfhcient. 
We expand the superconducting order parameter in 
layer n as ip n (r) = QY,q=o v n,q<P q ,N{0,4>) witn Q = 
($ k B T / pdoB) 1 / 4 and measure lengths in units of the 
magnetic length l m = (<&o/2ttB) 1 / 2 . In these units 
R={N/2f/H m . 

Within the LLL approximation, the Hamiltonian func- 
tional in Eq. (1) becomes a function of the complex co- 
efficients {v n ^ q } = v given by 
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where the effective temperature-field parameter is 20 



doQ 2 ( (rr , , \e\hnoH 
a T = - — — | ce(l ) 



k B T 



m± 



(4) 



[Note that ax = corresponds to the mean- field H C %(T) 
line and ax = — oo corresponds to T = 0.] The parame- 
ter that determines the coupling strength between adja- 
cent layers is then just 77 = H 2 /(2j 2 s 2 \eBh + m±a(T)\), 
where 7 = (m\\/m±) 1 / 2 is the anisotropy parameter. In 
the second term of Eq. (3) we have used 

N 
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xO(p-q)Q(N-p + q), (5) 

with 0(a:) the Heaviside step function, equal to zero for 
x < and one for x > 0. The Hamiltonian describing 
our model of a layered superconductor within the LLL 
approximation is that in Eq. (3). Note that W is a func- 
tion of the complex variables {v n , q } and depends on the 
temperature, Qr, and layer- coupling, rj, parameters. 

The equilibrium properties of this system are deter- 
mined by its partition function 

Z(T,H) = f Y[v^ n vr n -DA C xp(-H^ n ,A]/k B T), 

whose value is basically controlled by the order parameter 
configurations that minimize the Hamiltonian functional 



TL. In the LLL approximation, which neglects fluctu- 
ations in the magnetic induction and restricts the order 
parameter (in each layer) to the set of functions spanning 
the LLL, the partition function can be written as 



Z(T,H) = j Y[dv n , q dv* niq exp[-H(v)/k B T}. 
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Thermal averages of the quantities of interest X(T, H), 
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are calculated by means of MC simulations using the 
standard Metropolis algorithm 21 . We calculate thermal 
averages by the Importance Sampling Method 22 . 

The physical quantities that we have focused on to ex- 
amine in-layer vortex correlations are the structure factor 
(in each layer) defined as 23 
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where {x; } arc the positions of the N pancake vortices 
(on that layer) and we have chosen to parametrize the 
wave vector in polar coordinates as k = (k x ,k y ) = 
(k cos (j),k sin cj)), and the rotationally averaged structure 
function (also in each layer) given by 
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A Lorentzian fit to the first peak of this function is in 
agreement with an exponential decay of vortex density- 
correlations in real space 16 . The length scale governing 
that decay, £_d, is proportional to the inverse- width at 
half peak, S^ 1 , of the Lorentzian. 

Numerical measurements of these quantities have been 
made for runs up to 1.92 x 10 6 MC steps 24 with sys- 
tems as big as 18 layers containing 72 vortices per layer. 
We have studied the range of effective temperatures 
— 13 < cut < 2 (while cooling) for inter-layer couplings 
r) = 0.14, 0.5, 1, 2, 10, 100. One of the most notable fea- 
tures of a freezing transition, which is the appearance of 
Bragg peaks in the structure factor S(k), is observed at 
an effective temperature ax — —4 (for r\ = 0.14) and 
at higher temperatures for rj > 0.5. Fig. 1 shows this 
structure factor as temperature is lowered (see caption). 
Fig. 1(b) corresponds to the temperature ax = —4 at 
which the vortex system is essentially crystallized. We 
notice that at and below ax = —5 (Figs. 1(c) to 1(f)) 
the Bragg peaks appear exactly at the same positions. 
This is consistent with a legitimate crystalline vortex- 
phase for these low temperatures. Another signature 
of the vortex-liquid freezing is seen in the rapid growth 
of vortex density-correlations, below ax = —3 (for 
the same coupling 77 = 0.14). This is shown in Fig. 2, 
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FIG. 1: Structure factor, S(k), at layer M/2 of the clean lay- 
ered superconductor with interlayer coupling r\ — 0.14 for 
(N, M) = (72, 18) as the temperature is lowered: ar = 
-3,-4,-5,-6,-7,-8. [(a),(b),(c),(d),(e),(f) respectively.] 
The central maximum at k = has been removed for clar- 
ity. Note the appearance of Bragg peaks at a temperature 
between qt = —3 and —4. 



where the scaled density correlation length £d/Ris plot- 
ted. Fig. 2(a) shows the abrupt increase of this length 
scale as a function of our temperature parameter olt 
for all system sizes studied. Fig. 2(b) shows the same 
characteristic length-scale, £d/R, plotted (in log-linear 
scale) against |cvt| 2 - We observe that an exponential fit 
(straight-line there) to the data is appropriate at high 
temperatures (i.e. at low values of |o;t| 2 ), but breaks 
down in the low temperature regime where correlations 
grow even faster than that. The exponential behavior, 
£d ~ exp(const|ar| 2 ), has been predicted by Moore 25 
in connection with bulk anisotropic-superconductors in 
the continuum limit at their lower critical dimension. In 
that investigation, however, this rapid growth of cor- 
relations is not attributed to a thermodynamic phase 
transition but just to a crossover. Returning to Fig. 
2(b), we remark that the very last points mark the onset 
of in-layer vortex correlations growing even faster than 
£d ~ cxp(const | or | 2 ) [compare with Fig. 2(a) for lower 
temperatures.] We believe this can occur because a ther- 
modynamic first-order transition takes place at the tem- 
perature aT,m where £d departs from the straight line in 
Fig. 2(b). This takes place at a T , m = -3.25 for r) = 0.14. 
In fact, comparison with the work of Hu and MacDon- 
ald 12 gives excellent agreement for the melting (freezing) 
transition at these values. 

To study pancake-vortex correlations along the axis 
perpendicular to the layers (the c-axis), we measure nu- 
merically the two-point correlation function 
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FIG. 2: In-layer density correlations at layer M/2 of the 
clean layered-HTSC as the temperature is lowered, (a) Scaled 
vortex-density correlation- length, (,d/R, for different system 
sizes at rj = 0.14. The dashed curve is just a guide to the eye. 
(b) The same characteristic length-scale, £d/R, as a function 
of \ar\ 2 in log-linear scale for the system (N,M) — (72, 18) 
with exponential fit (solid line). [The arrow marks the tran- 
sition temperature for r\ = 0.14, ar,m = —3.25.] 
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FIG. 3: Angular vortex-positions in all layers of the pure 
layered-superconductor with interlayer coupling rj = 0.14 for 
(N, M) = (72, 18) as the temperature is lowered: olt = 
-3,-4,-5,-6,-7,-8. [(a),(b),(c),(d),(c),(f) respectively.] 
Equal gray-tone represent the same layer. The appearance of 
pancake- vortex domains signals the beginning of the coupled 
phase [which also appears at a temperature between ar = — 3 
and —4.] 



Here the overline denotes a spatial average over the 
surface of the sphere of radius R and the prefactor 
is just a "normalization" factor. We impose periodic 
boundary conditions on the order parameter: ipM+p = 
ip p , p — 1,2, ... ,M. Appearance of long-range inter- 
layer correlations can be observed directly from plots 
of the vortex positions in all layers as the temperature 
is lowered below the effective temperature a^m (sec 
Fig. 3) and are confirmed by the change from expo- 
nential to algebraic decay of the phase correlation func- 
tion, Re[Cp(m)], in Fig. 4(a). In the high-temperature 



4 




123456789 

m 



(a) 



6 8 10 12 

|oc r | 2 

(b) 



FIG. 4: Inter-layer phase correlations on the pure layered- 
HTSC as temperature is lowered, (a) Log-linear plot of the 
inter-layer phase-correlation-function, Re[Cp(m)], as a func- 
tion of the "distance" m along the c-axis. Exponential decay 
is observed above (open symbols) and about (crossed sym- 
bols) the transition temperature, ar,m, but power- law de- 
crease (filled symbols) below the transition. In each case 
the mentioned curve is fitted, (b) Scaled phase-correlation- 
length, £p/M, vs |o!t| 2 in log-linear scale for the system in (a) 
and also with exponential fit. [The arrow marks the coupling 
transition temperature, which to our resolution coincides with 

OtT,m-] 



regime, where the phase correlation function decays ex- 
ponentially as Re[Cp(m)] <~ cxp(— m/£p), we extract 
the phase correlation- length, £p, and plot it (scaled by 
M) against \ar\ 2 (Fig. 4(b)). We remark that, anal- 
ogously to £d, Ref. 25 also predicts the behavior £p ~ 
exp (const | evp | 2 ). This form gives a consistent fit to the 
data just above ar, m (solid line in Fig. 4(b)), but for 
olt < OLT, m the length scale £p grows faster than this 
exponential. We attribute the even more-rapid increase 
of £p for cvp < otT.m to an inter-layer coupling transi- 



tion which takes place also at ap.m- This criterion gives 
us good agreement also with numerical values of Ref. 15 
where the estimated coupling-transition temperature is 
about aT, m = —3.5 for 77 = 0.14. In that investigation, 
however, no freezing transition is observed! That is the 
fundamental difference argued in this paper. 

To summarize, we have investigated correlation func- 
tions for in-layer structure of the pancake- vortex system 
and for the coupling of the order parameter in different 
layers in the S-LD-LLL model. Results from MC sim- 
ulations, measuring the structure factor in every layer, 
suggest that there is an effective temperature aT, m below 
which the pancake-vortex system freezes to a crystalline 
phase but above which the system behaves as a liquid. 
This is signaled both by appearance of Bragg peaks in 
the vortex structure-factor S'(k) and the rapid increase 
of the in-layer vortex-density correlation length, £p . Si- 
multaneous to this freezing transition there appears a 
coupling of the order parameter in the different layers 
which is shown directly from measurements of the vortex 
positions in every layer and is confirmed by the existence 
of long-ranged correlations in the order parameter of the 
distant layers. Explicitly, the exponential decay in corre- 
lations along the c-axis above the transition temperature 
C(T,m changes to algebraic decrease in the low tempera- 
ture regime, cvp < ar,m- The transition is of first order 
and a.T,m depends on the interlayer coupling-parameter 
77. The LD-LLL model has shown 12 ' 15 consistency with 
the experimental melting-curve in YBCO. 
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